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Abstract 

We theoretically study electrically tunable magnetoplasmons in a monolayer of silicene or ger¬ 
manene. We derive the dynamical response function and take into account the effects of strong 
spin-orbit coupling (SOC) and of an external electric filed E z perpendicular to the plane of the 
buckled silicene/germanene. Employing the random-phase approximation we analyze the magne- 
toplasmon spectrum. The dispersion relation has the same form as in a two-dimensional electron 
gas with the cyclotron and plasma frequencies modified due to the SOC and the field E z . In the 
absence of SOC and E z , our results agree well with recent experiments on graphene. The predicted 
effects could be tested by experiments similar to those on graphene and would be useful for future 
spintronics and optoelectronic devices. 

PACS numbers: 71.45.Gm, 71.70.Di, 73.43.Lp, 78.30.Na 
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I. INTRODUCTION 


Since its realization as a truly two-dimensional (2D) material, graphene has attracted 
much interest, both due to fundamental science and technological importance in various fields 


.a. 


However, the realization of a tunable band gap, suitable for device fabrications, is still 
challenging and SOC is very weak in graphene. To overcome these limitations researchers 
have been increasingly studying similar materials. One such material, called silicene, is a 
monolayer honeycomb structure of silicon and has been predicted to be stable jsj]. Already 
several attempts have been made to synthesize it [l], A similar material is germanene. 

Despite controversy over whether silicene has been experimentally created or not 0 , it 
is expected to be an excellent candidate materials because it has a strong SOC and an elec¬ 
trically tunable band gap [Is- 8]. It’s a single layer of silicon atoms with a honeycomb lattice 
structure and compatible with silicon-based electronics that dominates the semiconductor 
industry. Silicene has Dirac cones similar to those of graphene and density functional cal¬ 
culations showed that the SOC gap induced in it is about 1.55 meV |], 7]. Moreover, very 
recent theoretical studies predict the stability of silicene on non metallic surfaces such as 


graphene 


0 . 


boron nitride or SiC 


10], and in graphene-silicene-graphene structures 11 ]. 


Besides the strong SOC, another salient feature of silicene is its buckled structure with the A 
and B sublattice planes separated by a vertical distance 2£ so that inversion symmetry can 
be broken by an external electric field resulting in a staggered potential js]. Accordingly, the 
energy gap in it and in germanene can be controlled electrically. Due to this unusual band 
structure, silicene and germanene are expected to show exotic properties such as quantum 


14 


m, 


spin- and valley-Hall effects js, 12, 13j, magneto-optical and electrical transport 
etc.. 

Plasmons are quantized charge excitations due to the Coulomb interaction and a very 
important aspect in condensed matter physics not only from a fundamental point of view 


but also from a technological one 


1 ( 3 - 20 ]. 


In the presence of a magnetic field they are called 


magnetoplasmons and have been extensively studied theoretically 2l] [25| and observed ex¬ 


perimentally 26-28] in graphene. The study of graphene (magneto)plasmons involves spatial 
confinement of light and enables them to operate at terahertz frequencies thus making it a 
promising material for optoelectronics. Next to graphene, which has a very weak SOC and 
no gap if not grown on a substrate, is silicene or germanene with strong SOC and a tunable 
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band gap. So far plasmons in them have been studied only in the absence of a magnetic 
field bd, 30], 31]. 

The purpose of this work is to study magnetoplasmons in silicene or germanene. We 
evaluate the dynamical nonlocal dielectric response function to obtain the magnetoplasmon 
spectrum within the random-phase approximation (RPA). In particular, we take into account 
the effect of strong SOC and of an external electric field E z applied perpendicular to its plane. 
Experiments can be done by incorporating the effects of SOC and E z similar to the recent 
ones 


26 


28] on gapless graphene. In Sec. II we present the basic formalism, in Sec. Ill 


the density-density correlation function, and in Sec. IV the magnetoplasmons. Results and 
their discussion follow in Sec. V and a summary in Sec. VI. 


II. MODEL FORMULATION 


We consider silicene or germanene in the ( x , y) plane in the presence of intrinsic SOC and 
of an external electric field E z applied along 2 axis in addition to a magnetic field B = Bz. 
Electrons in silicene obey the 2D Dirac-like Hamiltonian 7, 8] 


H v ’ s = v F (rja x U x + cr y Il y ) + r]s\a z + Va z . 


( 1 ) 


Here 7] = 1(—1) represents the K ( K') valley, V = 2 IE Z is the potential due to the uniform 
electric field E z , 21 = 0.046 nm is the distance between the two sublattice planes, and 
A = 4 meV the SOC. For germanene we have 21 = 0.066 nm and A = 43 meV. Also, 
cr x , a y , a z ) are the Pauli matrices that describe the sublattice pseudospin, v F the electron 
Fermi velocity, and s = +1(—1) the up (down) electron spin. Further, II = p — eA is the 
canonical momentum and A the vector potential that yields B = Bz] we use the Landau 
gauge A = (0,Bx, 0). After diagonalizing the Hamiltonian we obtain the eigenvalues 

K' = ± AV" + p]V K‘ = -iv e (2) 


where Vg = V + £A and £ = r/s. The corresponding eigenfunctions are 


e ik yy ( i(x) \ ^ +;S _ e^_ [ 0 | ( M*) | 

V^y y D n s (f) n (x) J ° \/%j y</> 0 (^) / ° VLy \ 0 / 


( 3 ) 


Here c u c = v F \/2eB/h, x = x — xo, xq = Pk y , l = \Jh/eB is the magnetic length, and 
L y the length of the silicene or germannene monolayer along the y direction. Moreover, 
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4> n (x) = e x2 / 2 H n (x)/^2 n n\y/nl and H n (x ) are the Hermite polynomials. C% ,s and D% s are 
the normalization constants 

cr = [(1 ± V ( /E?)/ 2] 1/2 , DT = [(1 =F V t IEr)/2] il2 . (4) 

The energy spectrum given in Eq. (2) is degenerate with respect to the wave vector k y . 
The eigenfunctions for the K' valley can be obtained from Eq. (3), by interchanging <f n and 
(j) n - 1 , and the corresponding eigenvalues from Eq. (2) with r) = —1. 


III. DENSITY-DENSITY CORRELATION FUNCTION 


i) Finite frequencies. The dynamic and static response properties of an electron system 
are embodied in the structure of the density-density correlation function which we evaluate 
in the RPA. The RPA treatment presented here is by its nature a high-density approximation 
that has been successfully employed in the study of collective excitations in 2D graphene-like 
systems both with and without an applied magnetic held 16]- 24]. It has been found that 
the RPA predictions of plasmon spectra are in excellent agreement with experimental results 


26]- 28]. Following this technique, one can express the dielectric function as 


e(q,uj) = 1-v c (q)U 0 (q,u), (5) 

where v c (q) = 2ne 2 /nq is the 2D Fourier transform of the Coulomb potential with wave 
vector q and n the effective background dielectric constant. The non-interacting density- 
density correlation function is obtained as 

n„(?,v) = i ;y [/(sr)-/w)]|<«'n- i,r ia>i 2 w 

n,n' ,k y ,ky 

x [£,r - K’ + + ilT\ 


where A is the area of the system and |a) = | n,rj,s,k y ). Here 7 is the the width of the 
energy levels due to scattering and is an infinitesimally small quantity in samples with high 


mobility 


28]. The matrix element in Eq. ( 6 ) is evaluated in the Appendix; the result is 


17' I e- 4 -' | o>|'- XAu) = & K ^,„{[CrCl,’}F n _^(u) + (?) 


where u = l 2 q 2 / 2. For n < n' we have [ F nn i(u )] 2 = ( n\/n'\)e U u n ' n [L^ n (u)} 2 and for 
n' < n the same expression with n and n' interchanged. The sum over k y in Eq. ( 6 ) can be 
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evaluated using the prescription (fc 0 = L x /2l 2 ) 


2tt 9s9v 

ky 


f k ° A 

I dky ~^~ 9 s 9 vt 

'-ko 


Dn 


( 8 ) 


where D 0 = 2i tI 2 , g s and g v are the spin and valley degeneracies, respectively. We use g s = 
g v — 1 in the present work due to the lifting of the spin and valley degeneracies in silicene 
or germanene. 

We now use the transformation k y —k y and the fact that E% s (k y ) is an even function 
of k y , see Eq. (2). Then if we interchange n and n! and perform the k y integration using 
Eqs. (7) and (8), we can write the non-interacting density-density correlation function as 


1 


n»( 9 ,w) = 


X 


«* - K‘ + fa + nr - (K‘ - K' + t*> + h) 


-1 


( 9 ) 


The real and imaginary parts of IIoA to) can be obtained from the identity l/{x ± i'y) = 
(p/x) =F i7r<5(x) where p denotes the principal value of l/x. The real part of Eq. (9) reads 

n 1 (q r , to) = -^-V' J nn > («)[/i(o;) + h (-w)], (10) 

A , 

n,n' 

with 

/.(•») = /(£?')/Kr - ey +m. (ii) 


while the imaginary part is written as 

n 2 (g,w) = J^2Jnn'(u)[l2(u) -/ 2 (-w)], 

with 


n,n' 


h(oj) = fwmhu + Er - ED. 


( 12 ) 


(13) 


Equations (10)-(13) will be the starting point of our treatment of magnetoplasmons. Their 
form makes clear their even and odd symmetry with respect to to. These functions are the 
essential ingredients for theoretical considerations of such diverse problems as high-frequency 
and steady-state transport, static and dynamic screening, and correlation phenomena. 

ii) Limit to — q = 0. The non-interacting density-density correlation function is obtained 
from Eq. (6) in the static and long wavelength limit, to = q = 0. Thus Eq. (6) becomes 


n 0 (o, 0) 


1 

A) 


E 

n,n',d= 


/(AT) - / W) 



(14) 
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where the summation over +/- represents electrons/holes. With the zero-temperature limit, 
this turns into a series of delta functions, 5(E F ± E% ,s ) 

6(E) = (T/n)/(E 2 + T 2 ), we arrive at 


32 


33]. Making the replacement 


n 0 (o, 0) — 


E 


(2 - V»)r 


2i,D « „sr± - B "*> 2 + r2 ’ 


(15) 


where T is the level width. Then the density-density correlation function is proportional to 
the density of states at the Fermi energy, fl 0 (0,0) = D(E F ). At finite temperatures though 


it is given by 
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33] 


n 0 (o, o) — 


r»+oo 


-df(E)/dE\ D(E)dE. 


(16) 


The density-density correlation function shows the lifting of the four-fold degeneracy at 
Ep = 0 (Dirac point) at zero temperature. At T = 0 and E F = 0 this function vanishes 
in the limit of zero SOC and E z , simply because it becomes the same as that of graphene 
at the Dirac point (E F = 0) with a completely filled valence band and completely empty 
conduction band. The corresponding carrier density vanishes and implies that no intrinsic 
graphene plasmons are possible (more generally, Dirac plasmons). This means that the 
screening is absent to linear order except for the renorrmalization of the dielectric constant 
term. However, when the Fermi level is away from E F = 0 or at nonzero temperature, the 
density-density correlation function shows doubly degenerate spin ad valley splitting of the 
Landau levels (LLs) and the linear screening is expected to become appropriate. Moreover, 
these results can be reduced to those for gappless graphene derived and discussed in Ref. 


33] (see Fig. 1) in the limit of zero SOC and E z . 


We show numerical results of Eq. (16) as a function of the Fermi energy in Fig. 1. We 
find that the n = 0 LL is split into four levels and all other LLs (n > 0) into two. The 
valley degeneracy is lifted by the application of the field E~ and the spin degeneracy by the 
SOC. This is consistent with the eigenvalues given by Eq. (2). We use B = 1 Tesla, T — 3 
K, and vary the field energy V = 2 IE Z and the SOC strength. The left panel is drawn for 
V = A = 0 (dashed curves) and V = 0 meV and A = 4 meV ( solid curves). The other two 
panels are for V = 7 meV and A = 4 meV; the middle panel is for the K valley and the the 
right one for the K' valley. 

In Fig. 1 the SOC and field split the LLs in two groups: in accordance with Eq. (2), 
rj, s = ±, we label them as + = +,+ = —, — and — = Every n ^ 0 LL is 
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FIG. 1. Static density-density correlation function in the long-wavelength limit q —>• 0 versus Fermi 
energy Ep. We vary the field E z and the SOC strength A. Left panel: the dashed and solid a 
curves are for V = A = 0 and V = 0 meV and A = 4 meV. The other two panels are for V = 7 
meV and A = 4 meV; the middle panel is for the K valley and the the right one for the K' valley. 
The degeneracy of the LLs is lifted. 


doubly degenerate in each group and consists of a spin-up state from one valley and a spin- 
down state from the other valley. The LL splitting between the two groups is symmetric in 
the valence and conduction band due to the symmetry in Eq. (2). The four-fold spin and 
valley degeneracy of the n — 0 LL is lifted by the SOC and electric held energy. 

iii) zero frequency. The static limit cu —» 0 of Eq. (6) is obtained with the help of Eqs. 
(7)-(8). In this limit Imll 0 (g,w) —> 0 and Eq. (6) gives 


n 0 M) 


1 v f(Ey) - f(E^) 
D 0 ^ EV - Ey 

n,n r 


Jnn' (^0 • 


(17) 


We show numerical results for IIo(g, 0) as a function of the wave vector q in Fig. 2. We use 
the parameters B = 5 Tesla, T = 10 K, and vary the held energy V = 2 IE Z and the SOC 
strength A. The black curves are for V = A = 0, the red ones for V = 0 meV and A = 4 
meV, and the blue curves for V — 10 meV and A = 4 meV. The solid and dashed curves 
pertain, respectively, to spin up and K valley and to spin down and K' valley. 

In the usual 2DEG the screening wave vector is independent of the carrier density but for 
graphene or silicene it is proportional to the square root of the density [if]]. First, in the limit 
of zero magnetic held B the static correlation function remains constant and equal to the 
electronic density of states up to the wave vector of q = 2kp] there are two contributions to 
it that stem from intraband and interband plasmons, respectively. In the large momentum 
transfer regime of Fig. 2, q ~ 5 (10 s m -1 , the static screening for the intraband case 


7 


















































































































decreases linearly with q, which is consistent with the case of gapless graphene in the limit 


of zero (see Fig. 2 of Ref. 


and Fig. 4 of Ref. 


32]) and finite |23J magnetic held. 


There is no possibility of zero-energy plasmon excitations in the intraband region (valence 
or conduction band). 



FIG. 2. Static density-density correlation function versus the wave vector q. We vary the electric 
field energy V and the SOC strength A. Black curves: V = A = 0; red curves V = 0 meV, A = 4 
meV; blue curves V = 1(1 meV and A = 4 meV. The solid and dashed curves pertain, respectively, 
to spin up and spin down in the K valley or spin down and spin up in the K' valley. Here we 
cannot distinguish between the black (or red) solid and dashed curves as there is no spin or valley 
splitting for the chosen parameters. 


We find a similar behaviour for finite B except in the small wave vector limit. In contrast 
with its behaviour at B = 0, the static correlation function tends to zero as —y 0, 0) oc q 2 


for finite B 


23]. This is due to the fact that the main contribution to it conies from the 


q = 0 excitations in the vicinity of Ep. Whereas at B = 0 there are q —> 0 excitations whose 
energy tends to zero, Ep now lies in a cyclotron gap between the highest occupied landau 
level 7i F and lowest unoccupied np + 1. This gap must be overcome by small-g excitations, 
such that its spectral weight approaches zero. The static correlation function also coincides 
with the density of states at Ep because the latter vanishes for finite fields B when Ep is 
in the gap. Further, the oscillatory behaviour of the static correlation function below 2 kp is 
due to intraband transitions, whether Ep is in the valence or conduction band (rip + 1, np). 





IV. MAGNETOPLASMONS 


Magnetoplasmons are readily furnished by the singularities of the function fli(< 7 , u), from 
the roots of the longitudinal magnetoplasmon dispersion relation obtained from Eq. (9) as 


1 - tUgJiMg.w) = 0, 


(18) 


along with the condition n 2 (q,uj) = 0 to ensure long-lived excitations 
is in excellent agreement with high-mobility graphene samples 28]. 


22 


23 


29 


30], which 


For weak damping the decay rate 7 , determined by Eqs. ( 10 ) and ( 12 ), is given by Eq. 
( 22 ) of Ref. jsoj]. Since we are primarily interested in the long-wavelength behavior of 
undamped magnetoplasmons, described by 7 oc n 2 (g, oj) = 0 , we treat them by solving Eq. 
( 18 ). With the help of Eq. ( 10 ) we End its roots are obtained by solving 

^2 


1 = 


kql 2 


’ S y^ J Jnn'{u) [h{u) + /i(—CU)] , 


(19) 


n,n' 


Using Eq. (11) we can write 

h M + h(~u) = 


2A V,S , 

n,n r 


h 2 ui 2 


(A- 0 2 


f(K’ s ), 


( 20 ) 


where A]// ; = E^f — E% s . Next we expand J nn '(u) to lowest order in its argument (low 
wave-number expansion). This amounts to considering only the n' = nil terms in Eq. (19). 
The inter-Landau level plasmon modes under consideration arise from neighbouring Landau 


levels, that is, from n' — n i 1. Then using the expansion 34] L l n (u) = Til 


for / > 0 and retaining only terms that are constant or linear in u we get 


m =0 


Jn,n-\- 1(^) —^ nuG c + {u + 1 )uG D , (21) 

Jn t n-i(u) -A (n — 1 )uG c i nuG°. ( 22 ) 

Here G c = (1 + ry)/ 2 , G D = (1 — ry)/ 2 , and 77 = V^/[H 2 uj 2 n + V 2 } 1 ! 2 . The factors G c and 
G D arise from the normalization of the eigenstates and in the limit A = E z = 0 are both 
equal to 1/2. 

To obtain the magnetoplasmon spectrum, we evaluate A^ s , for n' — n i 1 . We hnd 

A7„ ±1 = ± &; c /(2 [n + (U/fc c ) 2 ] 1/2 ). (23) 
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Substitution of Eqs (20)-(22) into Eq. (19) yields 


1 


e 2 q 

K 


E 


\K 


n,nd=l I 


nv - (a y„ ±1 )- 


f(A 


77,5 \ 

n ) 


(24) 


For inter-LL excitations near the Fermi energy Ep we can approximate n by hf in 
where rip is the LL index corresponding to Ep. This gives 


n 2 u 2 = (h 2 uj 2 c /(A[n F + (v^/huj c ) 2 ])) 
2e 2 q [rip + (V^/hw c ) 2 ] 1/2 


(25) 


x 


1 + 


K,h0J r 




With Ep in the conduction band (Ep = H 2 uj 2 np + V 2 ) Eq. (25) can be expressed as 


2 ~ 2 I ~ 2 

U — LU C + UJ , 


(26) 


where 

Cj c = u c [fru c / (2 (K 2 oj 2 c n F + V 2 ) 1 ' 2 )}, (27) 

and 

Cj p = ujp [vp/(h 2 uj 2 np + V 2 ) l/i ] , (28) 

1 /2 

with uj p = [e 2 q7Tn c /n\ ~ and n c = J2f(E% ,s )/( ti£ 2 ) the 2D carrier density. 

n 

It is interesting that Eq. (24) can be applied to the usual 2DEG for which A^ ±1 = ±Huj c . 

Then we obtain again Eq. (26) with Cj c and Cj p replaced, respectively, by u j c = eB/m and 
1 /2 

bJ v = [e 2 g7rn c /fi;] ", that is, the well-known plasmon dispersion relation. One can also take 



(26), namely the y /q behaviour, is common to 2D electron gas systems while the carrier 
density dependence is characteristic of the linear-in-fc dispersion relation of massless Dirac 


fermions, for which Ep = tiVy/ivn c . However, in the present case we can see the effects of 
gapped silicene or germanene with massive Dirac fermions and spin/valley splitting due to 
the combination of the SOC and the electric field E z . 


10 






















V. DISCUSSION OF RESULTS 


A closer analytical examination of Eq. (26) shows the following aspects of the gapped 
magnetoplasmon spectrum. If we set E z = 0 in Eq. (26) we obtain a SOC-induced, small- 
gap magnetoplasmon spectrum. Increasing E z , we obtain a larger gap, splitting and tuning 
of plasmons in silicene by combining it with the SOC. If we use a held E z comparable to 
the SOC strength A, then we expect splitting of the magnetoplasmon modes due to the 
combination of the two in the quantity V + r]s\. With further increase in E z , e.g., E z = 2A 
we can see an enhanced spin and valley splitting of the magnetoplasmon spectrum due to 
the V + r]s A factor in Eq. (21). Moreover, we note that the realization of topological phase 
transitions could also be observed in the magnetoplasmon spectrum if we take E z zero or less 
than A (spin-Hall regime), comparable to A (semi-metallic regime), and then twice A (valley- 
Hall regime). The spin-Hall regime is a topological insulator while the valley-Hall one is a 
rand insulator. For B —> 0 these transitions are consistent with recent plasmon predictions 


29|,l30]. Below we consider the effect of an external held B using the parameters 7. 26 30 : 
q = 7r/100 nnr 1 , vp = 0.5 x 10 6 m/s, A = 4 meV for silicene (43 meV for germanene) on 
SiC with dielectric constant k ~ 4 (different values do not qualitatively affect the results), 
and carrier density n c = 0.5 x 10 16 m -2 giving Ep = 41.3 meV. 

The changes in the density of states D(E) discussed in Sec. Ill and the approximations 
used to obtain the magnetoplasmons are rehected in the dependence of Ep , e.g., on the 
mangetic held. At finite temperatures the 2D carrier density n c is n c = J~D(E)f(E)dE , 
with D(E) for the LL spectrum obtained as 


DM = 2- H -£’’*) + £ i(B-sr)) 

0 V n=l,r],s ) 


(29) 


the factor 1/2 refers the fact the degeneracy of the zero LL is half that of the other LLs. 
Using Eq. (29) the result for n c becomes 


n, = 


Dn 


£ /«' 


rj, s 


7 i=l, 77, s 


f(ET 


(30) 


For hxed carrier density, this determines Ep implicitly by solving numerically Eq. (30). 
We show the resulting Ep, as a function of the magnetic held B in Fig. 3, for V = 10 meV, 
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FIG. 3. Fermi energy as a function of magnetic field for fixed values V = 10 meV, A = 4 meV, and 
n c = 0.5 x 10 16 m -2 . The temperature is varied such that T = 10 K (solid) and T = 5 K (dotted). 

A = 4 meV, T = 10 K, and n c = 0.5 x 10 16 m~ 2 . Ep remains constant for low B below 2T, 
that is, in the limit of large n; above this value we see the jumps as Ep crosses the LLs. 




).0 0.2 0.4 0.6 0.8 1.0 
B (Tesla) 



1.0 0.2 0.4 0.6 0.8 1.0 
B (Tesla) 


FIG. 4. Band structure of silicene as a function of the magnetic field B. The blue curves correspond 
to A = V = 0 and the black and red dotted ones to E~\ and E~, respectively. The black dotted 
line shows Ep vs B evaluated numerically using Eq. (30). In the left panel we cannot distinguish 
between the blue and red dotted curves because at A = V = 4 meV the gap is zero for the (V — A) 
curves. In the middle panel we see a clear degree of spin and valley splitting for A = 4 meV and 
V = 10 meV. The right panel, for A = 4 meV and V = 15 meV, shows a significant degree of spin 
and valley splitting by electrical tuning. 


We present the eigenvalues given by Eq. (2) as a function of the field B for fixed values of 
A and E~ in Fig. 4. We also include the Ep versus field B curve (dotted line) for comparison 
and further discussion. We find the following: (i) In the limit of A = £E Z = 0 (blue curves), 
we obtain the \/~B dependence of the LL energies. In contrast, for finite A and variable E z 
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FIG. 5. Magnetopasmons as a function of the magnetic field B for a fixed Ep = 41.3 meV. Blue 
curves correspond to A = V = 0.Black and red dotted curves represent and Huj-, respectively. 
In the left panel, with A = V = 4 meV, we cannot distinguish between the red and blue dotted 
curves because at equal amount of V and A energies, the gap is zero for the (V — A) curves. The 
middle panel, for A = 4 meV and V = 10 meV, shows a clear signature of spin and valley splitting. 
The right panel, for A = 4 meV and V = 15 meV, shows a significant spin and valley splitting by 
electrical tuning. The colour code is the same as in Fig. 4. 




FIG. 6. Magnetopasmons as a function of the magnetic field B for a fixed Ep = 26 meV. The 
parameters and curve marking are the same as in Fig. 5. 

(black and red dotted curves), the energies of the lower LLs grow linearly with B rather 
than with \/~B because of the massive Dirac fermions in silicene or germanene. (ii) The 
combination of the held energy V = 2 IE~ and A splits the LLs in two groups designated as 
E±, with E+ = E+ ,+ = E~~ and E~ = E~ ,+ = E+’~. (iii) The energies of the two groups 
of LLs in the valence or conduction band have not only different slopes versus B but also 
shift rigidly for B —>■ 0 due to the finite band gap either by A or by the held E~. However, 
every n ^ 0 LL is still doubly degenerate in each group, consisting of a spin-up state from 
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one valley and a spin-down state from the other valley. A crossing occurs between the two 
groups, which is symmetric in the valence and conduction band due to the symmetry in Eq. 
( 2 ). 

In Fig. 5 we show the magnetoplasmon spectrum as a function of the held B for fixed 
Ep = 41.3 meV. For comparison with graphene experiments 28], we show numerical results 
using Eq. (26) for A = V = 0 (blue curve). These results agree well with Eq. (1) and Fig. 


2 of Ref. 28], exhibiting dependence on VB, if we replace vp = 0.5 x 10 6 m/s by its value 
in graphene vp — lx 10 6 m/s. In the middle panel, for finite A and V = 2 IE Z , we found 
two curves, the red dotted (V — A) and black (V + A) showing a spin and valley splitting. 
The red dotted curve is the same as the blue one and we can’t distinguish between the two 
because the gap for the red dotted line vanishes due to A — 2 IE Z = 0. As the gap due to 
A and V is small and we are in a highly doped regime, we can see a split between the red 
dotted and black curve for two magnetoplamon modes fko± defined as fko + = huj + ^ + = tvjj_ _ 
and hhj_ = ho;_. + = Hui + ^. Increasing V = 2 IE Z = 10 meV, we see an enhanced splitting 
between red dotted and black curves for fixed A = 4 meV (middle panel). Here the blue 
and red dotted curves are weakly separated as the gap vanishes for the blue and red dotted 
curves V — A = 6 meV. With further increase in V, V — 15 meV, we obtain a further 
enhanced splitting between the black and red dotted curves of the magnetoplasmon modes 
as shown in the right panel for fixed A = 4 meV. We also note that the blue and red dotted 
lines are well separated as gap is zero for the blue and V — A = 11 meV for the red dotted 
curve. 

We contrast our results with those of recent graphene experiments on high-mobility or 
weakly doped samples 28], in the limit A = 2 tE z = 0, by further decreasing the Fermi 
energy close to the Dirac point. First, we show the magnetoplasmon spectrum as a function 
of the held B for Ep = 26 meV in Fig. 6. We found a clear splitting between the black and 
red dotted curves for A = 2 IE Z = 4 meV (left panel) as in the left panel of Fig. 5. As A 
and E z are small and we are in a weakly doped regime, we can see a strong splitting for the 
magnetoplamon modes Hu±. Again here we cannot distinguish between the red dotted and 
blue curves for the same reason as in Fig. 5. Upon increasing V, e.g. to V = 10 meV, we 
see a large splitting between the red dotted and black curves for fixed A = 4 meV (middle 
panel). We can weakly distinguish between the blue and red dotted curves here since the 
gap is V — A = 6 meV for the red dotted curve and zero for the blue one. With further 
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increase in V, V — 15 meV, we obtain a significant splitting between the red dotted and 
black curves of the magnetoplasmon modes as shown in the right panel. Here we also note 
that the blue and red dotted curves are well separated compared to those in the right panel 
of the Fig. 5. Again the results exhibit a square-root dependence on B and agree with recent 
graphene theory 
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25] and experiments 28] in the limit A = 2 £E Z = 0 provided we use 


v F = lx 10 6 m/s. 

The experimentally observed 281, \/~B dependence of the spectrum referred to above, in 
the limit A = V —> 0, applies to high-mobility weakly-doped, graphene samples, cf. Fig. 6. 


For highly-doped samples 


26|, 


27] though that involve values of E F A, V, with E F of the 


order of 200 — 300 meV, the magnetoplasmon gaps and spilttitings reported above will be 
very difficult to achieve as they would require unrealistically high values of V. Notice though 
that our analysis for silicene also holds for germanene, a monolayer of germanium, which 
has a much stronger SOC than silicene |7, 8], A ~ 43 meV. In both cases the predicted gaps 
and spilttings are sizeable for E F not too far from the Dirac point. 

Another feature of our results is the magnetoplasmon gaps. Although not yet experimen¬ 
tally confirmed, the SOC induced gap in silicene is about 1.55 meV [6,7] and is expected 
to be observed using existing experimental techniques. In the present work on electrically 
tunable magnetoplasmons in silicene, we have obtained a gap of about 1 meV in Fig. 5 and 
12 meV in Fig. 6 tuned by an external perpendicular electric field, which can be further 
enhanced by increasing this electric field and lowering the Fermi energy of the system close 
to the Dirac point. We believe that this gap can be observed in experiments similar to those 
on high-mobility graphene samples studying magnetoplasmons [28]. 

A possible extension of our work would be to include an in-plane electric field and study 


magneto-e 
in Ref. 


ectric-plasmons. One could then use the eigenfunctions and eigenvalues derived 


36] for E z = A = 0 as a starting point. 


VI. SUMMARY 

We showed electrically tunable effects in the magnetoplasmon spectrum of silicene and 
germanene due to the spin and valley polarization. Employing the RPA and including 
the effects of SOC and of an external electric field, we found a significant splitting of the 
magnetoplasmon spectrum. Our results agree well with graphene theory and experiments 
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in the limit of vanishing SOC and electric held provided Ep is not too far from the Dirac 


point, that is, for weakly-doped graphene samples 


28], if we use graphene’s value for vp. We 


expect that experimental studies of these novel phenomena in silicene, similar to those of 
Ref. [28], will be very appropriate since they directly bear on the many-body properties 
of silicene or germanene. Encouraging in this direction is the very recently reported local 
formation of high-buckled silicene nanosheets realized on a M 0 S 2 surface [35]. 

Electronic address: lm.tahir06@alumni.imperial.ac.uk 

Appendix A 

Below we outline the derivation of Eq. ( 8 ). The factor J a , a '{u) in Eq. (7) is given by 


J a ,a'iu) = (a'| |e- iq - r | |a ) 2 = (a'| Ie“ ici r I |a) x <a| |e iq ' r | \a '), 


(A.l) 


where |a) = |s, n, 77 , k y ). Using the eigenfunctions given by Eq. (3) Eq. (A.l) takes the form 


(a'\ |e _iq ' r | | a) = / dye^-^F / dx 

y n.rj.s J 

’ n —rvi 


-icy cj) n -i(x) 


D —iq x x 


-icy<j> n - i(x) 


DyM*) 

(A.2) 

where the superscript T denotes the transpose of the column vector. With the help of the 
identity (1 /L y ) f dy e l ( ky ^ ky ~ qy ) y = 5k' y ,k y -q y we can write Eq. (A.2) as 

(cr | |e q | |a) dk' y ,k y —q y ^ ^ E n / n ( q X: ky Qy, ky ) T Eh'— i,ri— 1 ( q x , ky q y , ky ) . (A.3) 


n,r],s 


Similarly, 


n,r],s 


Combining Eqs. (A.3) and (A.4), we arrive at 

Jn,n'(F) ^kyiky—q y ^ y ^n',n ( ky Qyt ky) X F n n ' (q X 7 kyi ky Qy) 


n,r],s 


(«| |c | \& ) &k' y ,ky—q y ^ ^ (yQxi ky , ky Qy') T F n _l n l _l [^Q X , kyi ky Qy ) . (A.4) 


T E n '— 1 yi, — 1 ( Qxi ky Qy , ky ) X F n —\ r i'_l (yQxi ky , ky Qy ) . (A.5) 

Now we proceed with the evaluation of F n t n (—q x , k y — q yi k y ). Using the explicit form of the 
harmonic oscillator functions 4> n (x) we have 
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oo 


D>h s nr/,si r „ 

TP . ( _ n b — n k ) — — - dX p -( x+l ( k v~<ly)) / 1 2 3 4 

”" ( s M MV* 2 " 2 "'n!n'! J 

—OO 

x H n , (X + l(ky - q y )) e~ iqxX e~( x+lky)2/2 H n (X + lk y ), (A.6) 

where X = x/l. Making the change Y = X + lk y + l {—q y + iq x ) /2 in Eq. (A.6) yields 

Z7 . ('_/r k —a k )— ^ n ' ^ r ~ u2 r il2 <lx(-<lv+ 2k v)/ 2 

- L n'n V l Lxi Hyi —: ° ° 

\/ Ti2 n 2 n n\n'\ 

OO 

x J dY e~ Y2 H n , (Y - l (q y + iq x ) /2) H n (Y - l (q y - iq x ) /2), (A.7) 

— OO 

where u = l 2 q 2 /2. The integral over Y is tabulated in Ref. 34, pp. 838 #7.377. The result 
for n < n' is 

F n > n {-q x , k y - q y , k y ) = (n!/n'!) 1/2 e -ul 2 ^ qx (- qy+ 2 k y)/ 2 


X 


l ( q y - iq x )/V 2 


u . 


(A.8) 


For n' < n, the result is given by Eq. (A.8) with n and n' interchanged. Using Eqs. (A.5) 
and (A.8) we arrive at Eq. (7), 

= |{a'| e-<^ |a>| 2 = D™\F m , («) + {CyC^F.-u,^ (u) }*, (A.9) 

with F nn i (■ u ) given after Eq. (7) in the text. 
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